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ABSTRACT 

We investigate the association between the accretion disk during radio plateau state and the follow- 
ing superluminal relativistic radio jets with peak intensity varies from 200 mJy to 1000 mJy observed 
over a period of five years and present the evidences of direct accretion disc-jet connection in micro- 
quasar GRS 1915-f 105 . We have analysed RXTE PCA/HEXTE X-ray data and have found that 
the accretion rate, niaccr, as inferred from the X-ray flux, is very high during the radio plateaux. We 
suggest that the accretion disk during the radio plateaux always associated with radiation-driven wind 
which is manifested in the form of enhanced absorption column density for X-ray and the depleted IR 
emission. We find that the wind density increases with the accretion disk luminosity during the radio 
plateaux. The wind density is similar to the density of the warm absorber proposed in extragalactic 
AGNs and Quasars. We suggest a simple model for the origin of superluminal relativistic jets. Finally, 
We discuss the implications of this work for galactic microquasars and the extragalactic AGNs and 
Quasars. 

Subject headings: accretion, accretion disk - X-rays -black hole physics - quasars -stars: individual 
(GRS 1915+105 ) 



1. INTRODUCTION 

Large superluminal outfiows/jcts which appear to 
move faster than light are widely studied in active 
galactic nuclei (AGN) & quasars, yet remain poorly 
understood (first AGN was discovered in 1943 by 
Carl Seyfert!). This is due to the long time scales 
associated with these systems as well as the central 
part of these systems is usually obscured by large 
amount of dust. Microquasars in our galaxy are stellar 
mass analogs of the massive black hole systems in 
quasars and AGNs. Microquasars are closer, smaller 
and show faster var ia bility that is easily observable. 
iMirabel fc Rodrigued l)1994j) discovered first micro- 
quasar GRS 1915+105 in our galaxy with superluminal 
jets. The black hole r nass in GRS 1915+105 is deter- 
mined to be 14+4 Mq l|Greiner. Cubv. fc: McCaughreaiJ 
12001). This microquasar shows exceptionally high 
variability in both X-rays and radio. Since 1996, 
rich X-ray variability of this s ource is o b served 
by RXTE (Morgan. Remillard. lg Greineil Il997t 
iMuno. Morgan, fc R.emillard. J,999) and by Indian 
X-ray Astronomy Experi ment, IX AE IjYadav et al.l 
11999(1 . iBelloni et al .' (2000) have classified the complex 
X-ray variability of GRS 1915+105 in 12 separate classes 
on the basis of their light curves and the color-color 
diagrams and suggested three basic states of this source 
namely low hard state, high soft state and low soft state. 

The radio emission from GRS 1915+105 can be 
broadly classified into three types: (a) steady radio jets 
(radio plateau state), (b) oscillations/baby jets (discrete 
jets) of 20-40 min duration in infrared (IR) and radio, 
and (c) large superluminal radio jets. The steady radio 
jets of 20 - 160 mJy flux density are associated with the 
canonical low hard X-ray state and observe d for e xtended 

m i200l . These 
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are optically thick compact jets of < 200 AU with ve- 
locity of 0.1-0.4 JDhawan. Mirabel, fc R,odriguezl200nl: 
iRibo. Dhawan. fc Mirabel l2004|) . The radio emission is 
correlated with the X-r ay emission as Ijradin ex L y^ fo r 
several different sources IjGallo. Fender, fc Poolev"2003"). 
^oolcv fc Fender ( 1997) observed radio oscillations with 
delayed emission at lower frequency. Simultaneous X- 
ray, IR and radio multi-wavelength observations pro- 
vided first major step in our understanding of disk- 
jet interaction and suggested that the spike in X-ray 
coincides with the beginning of IR flare and it has 
self absorbed syn chrotron emission assoc i ated with adi- 
abatic expansion l)Eikenberrv et alJll998t IMirabel et alJ 
Il998|) . These are also compact jets with velocity 
/3 ~ 1 rDha wan. Mirabel, fc R.odrigued l200?itl . It is 
further suggested that the spike in X-ray is associ- 
ated with the change of the X-ray state due to ma- 
jor ejection episode l|Yadavl 12001(1 . This is consis- 
tent with the suggestion that it is the coronal ma- 
terial and not the disk material (or any other accre- 
tion flow associated with the l ow hard state) which 
is ejected prior to ra d io flares ( Ran fc Greiner 200^ ; 
Vadawaie et abl 1200^ iFender. Belloni. fc Gaila ,2QQ4; 



Rothstein. Eikenberrv. fc Matthewai2005(l . 



The relativistic superluminal jets with up to 1 Jy 
flux density have steep radio spectrum and are ob- 
served at larg e distances few hundred AU to 5000 AU 
from the core IMirabel fc Rodriguedll994HFender et al. I 
IT999t iDhmvan. Mirabel, fc R,odriguedl2000n . These ra- 
dio jets are very energetic with luminosity close to the 
Eddington lumi nosity, L^rfrf and have been observed in 
several sou rces (Fender et al. ||1999: Hicll ming fc RupenI 
IT995t (wTet al...2002: .Orosz et al.. .20011 Progress in 
our understanding of these jets specially their con- 
nection to accretion disk has been slow. A semi- 
quantitati ve model is proposed for the disk-jet couplinj 
((Fender. Belloni. fc Ga,1]dl200l IFender fc Be]lonill200 
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Fig. 1.— GBI 2.25 GHz radio data for 1999 January 30 (start at MJD 51204.7) and 1997 October 30 (start at MJD 50750.6) large 
superluminal radio flares in the top panel. Start of the flares is offset to 0. It includes preplateau flares, plateau state and the large 
superluniinal radio flare. The dotted line connects the data points of 1999 January 30 superluminal radio flare to guide eyes. The radio 
spectral index is shown in the bottom panel. 



which although puts all available observational details 
in perspective but the scope to gain new insight or 
new inform ation about disk-jet interaction seems lim- 
ited IjRemillard 2005) . The physical connection between 
X-ray emission and the superluminal flares has been 
the hardest to understand. However a general trend 
is emerging that the plateau states are generally fol- 
lowed (and perhaps preceded) by superluminal flares 
IjFender fc Belloni 2004) . 

In this paper, we investigate the association of large 
superluminal jets with the radio plateaux. We have 
analysed the available RXTE PCA/HEXTE X-ray data 
during radio plateaux and the radio flare data from the 
Green Bank Interferometer (GBI). We describe our se- 
lection criterion and discuss in detail the morphology of 
superluminal radio flares. We present tight correlations 
between parameters of the accretion disk and the super- 
luminal jets. We suggest that the accretion disk during 
radio plateaux always associated with radiation-driven 
wind and corroborate it with supporting evidences. Fi- 
nally we provide a simple model for the origin of super- 
luminal jets and discuss its implications for galactic & 
extragalactic radio jet sources. 

2. OBSERVATIONS AND ANALYSIS 



All the radio flares observed in GRS 1915+105 can be 
broadly put into two groups on the basis of their flux, 
radio spectrum and spatial distribution; (1) the super- 
luminal flares (200-1000 mJy) which have steep radio 
spectra and are seen at large distances (> 240 AU), 
and (2) all other flares (5-360 mJy) which include the 
preplateau flares, radio oscillations & discrete flares and 
the steady radio flares during the plateaux. All these 
flares have flat radio spectra and are observed close to 
the compact object (< 200AU). It is the superluminal 
flares which are hard to understand. The other flares 
in the latter group are studied in detail and under- 
stood fairly well (this is discussed in the last section). 
In the internal shock model for superluminal radio jets 
in microquasars, an exponent ial decay is suggested once 
the s hock is fully developed IjKaiser. Sunvaev. fc Spruia 
120001) . It is also found that the radio pla teau is always 
associated with superluminal radio flare fFcnder et al. I 
.1999; Klei n- Woh ct al. 2002; Vadawale et al. 200ai) . We 
searched 2.25 GHz GBI radio monitoring data during the 
period from 1996 December to 2000 April and selected 
radio plateaux and the following radio flares which de- 
cay exponentially. The preplateau radio flares do not 
follow exponential decay and are like the oscillations and 
discrete flares but closely spaced in time as described 
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in the next paragraph. Seven flares are found to sat- 
isfy our criteria for superluminal flares and are listed in 
Table Q] One more radio flare on 2001 July 16 which 
was observed by the Very Large Baseline Array (VLBA) 
and Ratan radio telescope is also added in Table Q] 
VLBA observations clearly showed an ejecta well sepa- 
rated from the core IjDhawan et al.l l20031. This plus two 
more flares out of eight flares listed in Table [T] sho w mov- 
ing ejecta well separated from the c ore ( Fender et al. I 
[19991 iDhawan. Mirabel, fc Rodrigue^ 2000). In Table U 
rise time is the total time taken by a radio flare to rise 
from plateau flux to the peak flux. The gaps in the 
GBI data put upper limit on the rise time (< 1 day). 
For 340 mJy flare on 1999 January 30 the rise time is 
0.25 day while for 710 mJy flare on 1998 June 3 the rise 
time is 0.3 day. These values are consistent with the re- 
ported values of 0.25 - 0.5 day f or different microquasars 
ijFender. Belloni. fc Gallcl l2004|) . The decay time con- 
stant is calculated by fitting exponential decay profile to 
the radio flare data. The contribution due to continuous 
radio plateau (like flares on 1998 April 13 & 30) or any 
other minor flare if any is removed prior to profile fitting. 
The decay time constant varies from 1.12 to 4.02 days. 
The superluminal radio fiare profile can be described as 
fast rise and exponential decay (FRED). Here fast rise 
(0.2-0.5 day) is in comparison to the decay time. For 
comparison the oscillations/discrete radio jets have rise 
and decay time in the range of 0.1 - 0.3 hr dMirab el et al.l 
119981 llshwara-Chandra. Yadav. fc Pramesh Rad l20o3r 
The integrated radio flux is calculated by integrating the 
fitted exponential function over duration three times the 
decay time constant. 

The typical sequence of events for a superluminal ra- 
dio flare is shown in Figure H for 550 mJy radio flare 
on 1997 October 30 (gap in GBI data at peak) and 
for 340 mJy radio flare on 1999 January 30. Start 
of the radio flares is offset to zero. A superlumi- 
nal event starts with small preplateau flares followed 
by steady long plateau, followed by superluminal radio 
flares. We have described the radio plateau state and 
the superluminal flares in previous section. The pre- 
plateau flares are studied in detail using GMRT and 
RYLE radio telescope dat a for 2001 July 16 super- 
luminal flare jl^hwa ra-Cha ndra. Yadav. fc Pramesh Raol 
120021 lYadav et al.ll2003^) . These are discrete ejections of 
adiabatically expanding synchrotron clouds with flat ra- 
dio emission (delayed emission at lower frequency). They 
are like the oscillations/discrete jets but closely spaced 
in time and hence produce overlapped radio flares. Ra- 
dio emission at 1.4 GHz and 15 GHz supports flat radio 
spectrum. 

We study X-ray properties during radio plateaux 
within the preceding week from the start of superluminal 
flares to avoid changes in accretion disk over the long du- 
rations of plateaux (from -7 to -2 in Figure ^). We also 
avoid last day of plateaux as radio data suggest rapid 
changes in the accretion disk. For the timing analysis, 
we used single-bit-mode RXTE/PCA data in the energy 
ranges 3.6 - 5.7 and 5.7 - 14.8 keV. Power density spec- 
tra of 256 bin light curves are generated and co-added for 
every 16 s. The power density spectra are normalised. 
We fit a model of power law -t- Lorentzian and have cal- 
culated quasi- periodic oscillation (QPO) frequency. We 
show in Figure|21the power density spectra in 0.09 - 9 Hz 
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Fig. 2. — Normalised power density spectra in 0.09 — 9 Hz fre- 
quency range observed on 1998 April 28 and 2001 July 11. Strong 
QPOs are seen at 1.4 Hz and 2.4 Hz in PDS observed on 1998 
April 28 and 2001 July 11 respectively. The first harmonics are 
also visible in spectra. 



frequency range observed on 1998 April 28 and 2001 July 
11. Strong QPOs are seen at 1.4 Hz and 2.4 Hz in PDS 
observed on 1998 April 28 and 2001 July 11 respectively. 
The first harmonics are also visible in the spectra. For 
spectral analysis, we used PGA standard-2 data (3 -35 
keV) and cluster data from HEXTE (20 - 150 keV). 
We added 0.5% systematic error in PGA data. Standard 
procedures for data reduction, response matrix and back- 
ground e stimation are followed usi ng HEASOFT (ver 
5.2). M uno. Morgan, fc RemiUardl |1999) have shown 
that a spectral fit using sum of multicolor disk blackbody, 
power law and a gaussian line during plateaux results 
high values of the inner disk temperature (> 3 kc V) and 
unrealistic accretion disk radius (< 2 km). Fuch s et alJ 
( 2003) also drew similar conclusion while analysing X-ray 
properties d uring a platea u state using INTEGRAL and 
RXTE data. iRau fc Gre iner ( 2003|) analysed large sam- 
ple of RXTE observations belonging to radio plateaux 
using multicolor disk-blackbody -I- power law spectrum 
reflected from the accretion disk. This model adequately 
describes the X-ray spectra of all the observations, but 
it yields extremely high values of reflection parameter 
inconsistent with other measurements. 

A spectral model of simple power law with multicolor 
disk-blackbody and a gaussian provides reasonable val- 
ues of temperature and inner radius of the disk if the 
energy spectrum is disk dominate d as is the case during 
class /3 IjMigliari fc Belloiii1l2003tl . During the plateaux, 
the source is in very high luminosity state (VHS) with 
power index F > 2 . The spectrum is dominated by 
the Gompton scattered emission (> 85%) rather than 
the disk. The comptonised emission has a complex spec- 
trum shape showing features from both thermal and non- 
thermal electrons (Gliris & Kubota 2005). It suggests 
the presence of an optically thick inner disk corona which 
drains energy from the disk and reduces the disk temper- 
ature. The presence of QPO in 1.4 - 2.6 Hz frequency 
range is consistent with the large corona (see Table ^. 
This will require additional component in the spectral 
model during the plateaux. A three component model 
(multicolor disk-blackbody -I- power law -I- a comptonised 
component (GompTT)) with a gaussian line at 6.4 keV 
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Fig. 3. — Correlation between the total 3-150 keV X-ray flux 
from spectral analysis of RXTE PCA/HEXTE X-ray data during 
preflare plateau state and the integrated flux of the superluminal 
radio flare. The solid line is the linear fit to the data points. 
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Fig. 4. — Correlation between the quasi-periodic oscillation 
(QPO) frequency from timing analysis of RXTE/PCA X-ray data 
during preflare plateau state and the decay time constant of the 
superluminal radio flare. The dotted line is the linear flt to the 
data. 



is reported to provide more realistic va lues of the inner 
disk temperature and the disk radius IjRao et alJ 1200(1 
iVadawale et al.lEoO,'^ . We have used this model for our 
X-ray spectral analysis during plateaux listed in Tablc^] 
The absorption column density, N/f is kept as free pa- 
rameter. The integrated X-ray flux for energy range 3 - 
150 keV and the derived values of Njy are listed in Ta- 
ble Q] The other best-fit X-ray spectral parameters are 
given in Table |21 along with the flux of the individual 
spectral components. The x^ of our spectral fits falls in 
the range 0.81 to 1.4 (for dof 74-94). The spectrum is 
dominated by the Compton scattering emission and the 
contribution from the disk is limited to less than 13%. 

3. RESULTS AND DISCUSSION 

It is widely believed that the coronal material (and 
not the disk material) is ejected prior to radio flares 
as discussed earlier. It is consistent with the obser- 
vation that radio flares are observed during transition 
from the low hard state to the high s oft state but never 
obser ved during opposite transition l)Klein-Wolt et alH 
|2002|) . All the radio flares observed close to the core 



in GRS 1915+105 can be further subdivided into two 
groups on the basis of associated change observed in 
the accretion disc; (1) the persistent radio flares dur- 
ing the plateaux with steady accretion disk (no change 
associated with radio flares) and, (2) all other radio 
flares which accompany clear change in the accretion 
disk. In the former case, inflow and outflow in the ac- 
cretion disk are in equilibrium. Both the accretion disk 
and the radio flare are in steady state and radio emis- 
sion is tightly correlated with the X-ray luminosity as 
discussed above. In latter case, the hard X-ray com- 
ponent (Compton component) is either partly or com- 
pletely affec ted. The large changes will resul t state 
transition llRothstein. Eikenberry. fc Matthews! I2005t 
iRau fc Greined 120031: lYadavll200H) . The accretion disk 
may come back to its initial state after some time or 
may remain in the new state for a longer duration. 
If the accreti on disk is in some periodic state like in 
class (3 or 9 JBelloni e t al J 120 001). it produ ces IR and 
radio oscillations. lEikenberrv et al.l l)1998t) suggested 
that the spike in X-ray during the class /3 is related 
to mass ejection and the ejected matter decouples from 
the accretion disk. It is further shown that the inten- 
sity of these IR and radio jets increases with the X- 
ray hardness ratio during the low hard state preced- 
ing the spike (Yadav 2001). For these oscillations and 
the preplateau radio flares, all the available experimen- 
tal data of time delay in emission at lower frequen- 
cies can be adequately explained using adiabatically ex- 
panding synchrotron emitting cloud ejected from the ac- 
cretion disk (Ishwara-Chandra, Yadav, & Pramcsh Rag 
|2002() . These results suggests that the ejected matter be- 
haves like an adiabatically expanding synchrotron cloud 
w hich has bee n decou pled from the accretion disk. 

iMuno et alJ 1)200 ll) have shown that radio emission 
during large superluminal flares is decoupled from the 
accretion disk. This view is further strengthened by 
the fact that superluminal flares appear at distances 
few hundred AU from the core JFender et al. 1 119991: 
iDhawan. Mirabel, fc Rodriguez! !200fl() . !Fender et all 
l) 19991) have observed superluminal flares after a radio 
plateau state. To investigate the association of superlu- 
minal flares with plateaux, we plot in Figure |31 the to- 
tal 3-150 keV X-ray flux during the preceding plateau 
vs the integrated radio flux of the superluminal flare. 
We drive a correlation coefficient of 0.99 suggesting a 
strong connection between the total X-ray flux during 
the preceding plateau state and the integrated radio flux 
of the superluminal flare. We also find strong correla- 
tion between ASM count rate vs plateau radio fiux and 
A SM count rate vs the in tegrated radio flux as discussed 
bv !Vadawale et al] f 2003) (not shown here). In Figure^ 
we plot QPO frequency from our timing analysis of X- 
ray data during the plateau state as a function of decay 
time constant of the following superluminal radio flare. 
It again shows a tight correlation with correlation coef- 
ficient of 0.98. The remarkable feature of our findings 
here is that the parameters calculated using completely 
independent spectral and timing analysis bring out clear 
connection between the accretion disk during the plateau 
state and the following superluminal radio flare. 

In the last column of Table Q we give calcu- 
lated value of Nh which ranges from 10x10^^ 
to 15x10^^ cm^^. These values are higher than 
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commonly used N/f ^^ 5 - 6x10^^ cm^^ for spec- 
tral analysis of GRS 19 15+105 by several authors 
llBelloni et al. Il997t 'Mun o. Morgan, fc RemillardI 
119991 lYada"^ 120011 feau & Greineil l20oHr RaT 

dio emission from atomic hydrogen & molecular 
hydrogen together puts line-of-sight column den- 
sity > 3x10^^^ cm"^ (iDickey fc LockmanI Il990t 
iDame Hartmann. fc ThaddeusI l200l|r iBelloni et aD 
(j2000() found for the first time eyidences for a variable 
Nh during the plateaux (class xi ^ Xs) ^^'^ suggested 
presence of a large intrinsic absorber. The value of N/f 
estimated to be 6x10^^ cm~^ for plateau state while 
for non -plateau low hard state, it is found to be 2x10^^ 
cm-2. iLee et al. I ||2002D have analysed CHANDRA & 
RXTE X-ray data during a radio plateau state with GBI 
radio flux of 20 mJy at 2.25 GHz on 2000 April 24. The 
total 2-25 keV X-ray flux is 1.89x10"* ergs cm"^ s"^ 
These plateau conditions are identical (in radio & X-ray 
flux) to the one preceding the superluminal flare on 2001 
July 16 listed in Tabled The K absorption edges from 
Fe, Si and others are detected in CHANDRA data which 
suggests the presence of a warm absorber in the vicinity 
of accretion disk. The N/f derived from Fe absorption 
edge is found to be 9.3tit'3 xlO^'^ cm"-^ for solar 
abundances which is in agreement with our calculated 
value of N// ~ 10.6;to'4 x 10^^ cm~^ for the superluminal 
flare on 2001 July 16 (see Figure |3J). The lower limit 
of bolometric luminosity L^oi ~ Lx = 6.4x10^* ergs 
s^^ which is 0.35 of Eddingto n Luminosity, L ^rfrf for 
a black hole of mass 14 M© IjLee et al. II2002|) . The 
corresponding lower limit of bolometric luminosity for 
the X-ray flux listed in Table [H falls in the range of 
0.35 - 0.48 of LEdd- ILee et al. I (,2002^) have suggested 
wind from the accretion disk as the source of observed 
enhanced N//. One can define wind (spherical) mass 
outflow rate, m^ind ^ 47rr^py(f7/47r) ^ 9.5xl0^*(il/47r) 
g s~^ where r is the radius (10 cm), p is the density 
of the absorbing material and v is the wind velocity 
(100 km s~^). The accretion rate, rhaccr = Lboi/(?yc^) > 
Lx/irjC^) ~ 7.1x10^* g s~^ where the efficiency r/ ~ 0.1. 
It shows that as covering fraction (ri/47r) approaches 
unity, the wind outflow rate, rriwind is comparable to 
the iTiaccr- During the low hard state, the efficiency of 
compact persistent radio jets is estimated to be ~ 5% 
llFendedl200lt iFender. Belloni. fc Gallcll200l . In this 
case, losses due to non-radiative process (like adiabatic 
expansion) are likely to dominate. Our estimate of 
fhwind above suggests that the major part of accretion 
energy goes to the wind during the plateaux. When 
GRS 1915-1-105 is accreting near Lsdd, the presence of a 
radiation-driven wind is always expected and the wind 
density should be a strong function of the disk luminos- 
ity. Our derived values of Nh show strong dependence 
on observed total X-ray flux with correlation coefhcient 
of 0.995 which is shown in Figure [S] The dotted line 
is a linear fit to our data (plus sign with error bars). 
The open circle with error bar shows the value derived 
from CHANDRA and RXTE data as discussed above 
which is in ag r eemen t with our data within error bars. 
iKotani et al. I l)2000l) have estimated Nh ~ lO^'' from 
ASCA data. Our estimates of Nh falls between these 
two values. 
The best-fit X-ray spectral parameters are given in Ta- 
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Fig. 5. — Correlation between the total 3-150 keV X-ray flux and 
the absorption column density, N^f (plus sign with error bars). The 
dotted line is the linear flt to our data. The open circle (with error 
bar) is the derived N^f for the solar abundances from CHANDRA 
and RXTE X-ray data by Lee et al. (2002) during a plateau state 
with 20 mJy radio flux which is similar to the one observed during 
the plateau preceding the superluminal flare on 2001 Julyl6 (for 
detail see text). 



ble 13 along with the flux of individual spectral compo- 
nents. The spectrum is dominated by the Compton scat- 
tering emission and the contribution from the disk is lim- 
ited to < 13%. The inner-disk temperature (kT.^) and 
the plasma temperature (kTg) are ^0.5 and 4-7 k eV re- 
spect ively which are in acceptable range IjVadawal e et alJ 
12003(1 . The derived inner-disk radii (Rm) are in the phys- 
ically plausible range of (10-25) Rg for 14±4 M© black 
hole. As the total X-ray flux rises from 1.75x10^* ergs 
cm~^ s"'^ to 2.42xl0~® ergs cm~^ s~^, the major part of 
it (close to 80%) goes to the COMPTT component which 
shows strong correlation with Nh (this is analogous to 
the solar wind originating from the dense solar corona). 
Our s pectral model d o es no t overestimate the value of 
Nh- iVadawale et al.l l)2003l) have analysed preplateau 
data of 2001 June 30 when source was in the low hard 
state (non-plateau, plateau started on 2001 July 03) and 
estimated N// ~ 1. 65 x 10^^ cm~^ whi ch is close to the 
value estimated by IBelloni et al.l l)200flD for non-plateau 
low hard state. We have also studied the X-ray class P in 
GRS 1915+105 using this spectral mode l which result s 
Nh ~ 7x10^2 cm-2 for xl close to one (|Yadav 1 12006). 
Similar N/f has been used earlier for spectr al studies of 
the X-ray class /3 l|Mighari fc Belloni ll200^ . 

Our calculated power density spectra shown in Fig- 
ure [3 also independently lend support to the presence of 
wind during the plateau state. In the high soft state, the 
PDS is a featureless power law of dP/dv oc v~°^ where 
a is ~ 4/3. On the other hand, the low hard state is 
dominated by the Compton scattering photon emission 
and the PDS develops bandwidth limite d noise, a flat- 
shoulder with a Q PO at 1-4 Hz ( Rao. Yadav. fc Paul I 
l200nHMeieni200,'iD . As the power of corona (Compton 
scattered emission) increases, the power in the PDS in- 
creases at higher frequency (i^ > 0.1 Hz) and the QPO 
frequency decreases. Figure [3 shows PDS spectra ob- 
served on 2001 July 11 and 1998 April 28. The Compton 
scattered flux increases from 1.5x10^* ergs cm^^ s^^ on 
2001 July 11 to 1.9x10-® ergs cm'^ s'^ on 1998 April 
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Fig. 6. — Correlation between the integrated radio fiux of the 
superluniinal radio flares and the absorption column density, Nj:/ 
calculated from X-ray spectral analysis of RXTE PCA/HEXTE 
data during preflare plateau state. The dotted line is a linear fit 
to the data. 



28 while observed QPO frequency decreases from 2.4 Hz 
on 2001 July 11 to 1.4 Hz on 1998 April 28. It is clear 
from Figure |21 that power at f > 0.2 Hz is less in the 
PDS observed on 1998 April 28 than that observed on 
2001 July 11. The fast variability is suppressed by pho- 
ton scattering in the enhanced wind on 1998 April 28 and 
hence reducing the power in the PDS at higher frequen- 
cies. rShaposhniko v fc Titarc huk (2006) have discussed 
the decrease in the PDS power observed in Cyg X-1 at 
higher frequency (i^ > 0.1 Hz) as wind increases. In Cyg 
X-3, the PDS power in the low hard state drops to be- 
low 10"'^ (rms/mean)^ Hz~^ at frequencies z/ > 0.1 Hz 
as dense wind fr om the comp anion always envelops the 
compact object (lYadav I l2006() . 

Kotani et al. ( 2000) analysed ASCA data of two 
GRS 1915+105 observations in the low hard state in 
1994 & 1995 and detected characteristic absorption lines 
which suggest large values of Nh- They have dis- 
cussed radiation-driven wind from the accretion disk as 
the source of enhanced Nh and estimated wind termi- 
nal velocity of the order of 10^ cm s""'^ for the val- 
ues of r and hx discussed above. They have specu- 
lated wind's role in the formation of superluniinal radio 
jets. The spectral analysis of energy spectra during the 
extended low soft state requires large Nh which lends 
further support for the radiation-driven wind from the 
accretion disk (Vadawalc. Rao. & Naik 2002). The evi- 
dences for a Seyfert-like warm absorber have been found 
in other microquasars like GRO J1655-40, GX 339-4 



and XTE J1650-500 beside GRS 1915-H05 ilMiller et al. 



2m 



D ITgM iLee et al. II200I iKotani et al 



20001) and superluminal flares have been observed in all 
these sources except XTE J1650-500. In Figure we 
plot the integrated radio flux of superluminal flare as 
a function of Nh calculated from the X-ray data during 
the preceding plateau state. It brings out clearly the con- 
nection between the superluminal flare and the enhanced 
Nh with correlation coefficient of 0.90. The dotted line 
is a linear fit to the data. The wind density for typical 
one AU sphere comes out to be ^ 10^ cm~^ which is 



similar to the value proposed in quasars l)Elvisll2000(l . 

The presence of a dense absorber close to the 
vicinity of the accretion disk should affect the IR 
emission during the radio plateaux. The source 
size of the steady compact jet during the ra- 
dio plateau state varies as a function of frequency 
which is consistent with self-absorbed steady outflo\y 
llDhawan. Mirabel, fc Rodr i(^zl2000t iFuchs et a'ni200l 
iFender. Belloni. fc G allo 20^^ It is similar to the 
self-absorbed discrete outflows in the case of oscil- 
lations/baby jets which produce delayed emission at 
lower frequencies llMirabel et al.lll998l: lEikenberrv et alJ 



T99Sl ITshwara-Ghandra. Yadav. fc Pra,mesh Ra,d I2002t 



ya: 
fl 



Ueda et aLll2002|) . In both the cases, the radio emis- 
sion is flat or inverted. The IR emission during radio 
plate a ux in GRS 1915+ 105 (Fuchs. Mirabel, fc GlareB 
120031: iFuchs et all l2003|) is closer to the IR emission 
in Cyg X-3 where wind fro m the companion always 
envelops the compact object J 9 s\^^- Burnell. fc Fended 
l200illKoch-Mira,mnnd et a,1. Il2002f). 

Our results support the inter nal shock model for the 
origin of superluminal flares ( K aiser. Sunvaev. fc Sprtira 
I2OOOI) . The internal shock should form in the previously 
generated slowl y moving wind from the accretion disk 
with /3 < 0.01 (|Kotani e t al. "2000') as the fast moving 
discrete jet with (3^1 (Dhawan. Mirabel. & Rodrigu^ 
I2000J) catches up and interacts with it. Both the com- 
ponents; slow moving wind and fast moving jet are re- 
lated to the accretion disk during plateau state and the 
strength fc speed of these two components will deter- 
mine the power of the internal shock. This aspect has 
been brought out clearly by the correlations presented in 
Figures Ol E fc El In our case Lorentz factors r2-ri ^ 
F-;. for these two component s l)Fender. Belloni. fc Gallol 
120041 iRees fc Meszar'o^ll994j) . The wind deposits large 
amount of energy as rhwind approaches rhaccr prior to the 
switching on of superluminal flare. GRS 1915+105 pro- 
duces upto 1 Jy superluminal flares with apparent veloc- 
ity, Papp in a range 1.2-1.7. It is a efficient internal shock 
scenario where not just velocity but total jet power has 
significantly increased (Fender. Belloni. fc Gallo 2004). 
This view is further strengthened by the fact that Cyg 
X-3 can produce superluminal flares with flux density 
upto 15 Jy with Papp ~ 0.69 where the companion star 
can provide denser wind. The internal shock model can 
easily ac commodate high jet power requirement > 10^* 
ergs s~^ l)Fender et al. II1999I) and can explain the shift- 
ing fro m thick to thin radio emission du ring superluminal 
flares (|Kaiser. Sunvaev. fc SpruiJl2000() . 

Radio plateaux occur when rhaccr, as inferred from the 
X-ray flux, is very high and hboi approaches l^Edd- As 
fTT-accr riscs, the accretion disk passes through an insta- 
bility zone (prior to plateau state) and produces pre- 
plateau flares (see Figure ^). Thereafter, it enters in 
steady plateau state. At the end of the plateau state, 
changes in rhaccr again produce instabilities in the accre- 
tion disk which produce post-plateau flare (flare which 
produce superluminal flare) and the source comes out of 
plateau state. It may be noted here that a post-plateau 
flare may not always end plateau state but it may change 
the level of the plateau as is the case of superluminal 
flares on 1998 April 13 fc 30 (plateau radio flux 91 mJy 
(highest flux in Table ^). The instabilities in the accre- 
tion disk may be triggered by either a change in rhaccr 
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or by some a djustment among different accretion flows 
IjYadav et alJ ll999). This view is supported by the find- 
ing that certain X-ray variability classes are observed 
preferably b efore and after the plateaux (low hard state) 
i|Naik et alJl20n2D . 

Our results presented here suggest only three types 
of radio emission in GRS 1915-1-105 ; (1) radio emission 
during the steady plateau state with /3 in a range 0.1-0.4 
(steady X-ray properties), (2) oscillations or discrete 
baby jets with (3^1 (state transition in X-rays), and 
(3) faint flares of 2-3 min durations with likely low (3 
hard X-ray spectrum changes but no st ate transition 
Rpthstein. Eikenberrv. fc Matthewll2005^ ). The super- 
luminal flares are the consequences of class 1 & class 
2 radio emissions when occur in this order (but not 
in the reverse order). Our description of superluminal 
flares here can explain why axis of superluminal flares 
differs by < 12° from the axis of AU s cale compact jets 
IjDhawan. Mirabel, k Rodriguez! 12000(1 . The wind is 
radiation-driven and is expected to originate from the 
inner part of the accretion disk while the compact jets 
are most likely to originate from the outer corona. The 
phase lag between the hard and soft X-ray photons in 
GRS 1915-fl05 changes sign durin g the plateau stat e 
and produce complex behavior l|Muno et alJ 1200 If) , 
^arnierc (2005) has suggested that anything located 
inside the inner radius of the accretion disk (in our case 
wind) that can absorb a small part of soft X-ray flux can 
explain the changing of sign of the lag and its complex 
behavior during the plateau state. Other consequences 
of our results is that one would not expect oscillations/ 
baby jet type of IR fc radio emission in Cyg X-3 
(no t observed so far IjDhawan. Mirabel, fc Rodrigued 
|2£)00)) as has been observed in 

GRS 1915+105 ("Mirabel et al." Iggs": 'Eikcnbcrrv et al." 
[1998; Ishwara-Chandra. Yadav. & Pramcsh Rao 2002). 
Any discrete jet (including oscillations) shoifld produce 
internal shock as wind from the companion is always 
present in the case of Cyg X-3 and therefore, such jet 
would always end up producing a superluminal jet with 
thin radio emission. 

The physics in all the systems dominated by black 
holes is essentially the same and therefore, it is al- 
ways tempting to compare microquasars with mas- 
sive extragalactic AGNs and quasars. The cor- 
relation between radio and X-ray luminosity dur- 
ing the low hard state ((Gallo. Fender, fc PoolevI I2003D 
can be extended t o the AGNs by including a black 
hole mass term llMerl oni^ Hei nz, fc di Matted 120031: 
iFalcke. Koerding. fc Markoffl I2004D Therefore the 

steady jets in microquasars during radio plateaux are di- 
rectly comparable with those in AGNs. iMarscher et al.l 
l|2002) have claimed to observe oscillation radio jet (~ 
1/year) in AGN 3C120. They related these radio jets 
to superluminal ejection events which can be under- 
stood in the frame work of the model discussed here 
(if we assume the presence of warm absorber). The 
results discussed in this paper suggest the propaga- 
tion of shock wave to produce the superluminal jets 
in microquasars, a interpretation usually favoured for 
jets in extragalactic systems. It is shown recently that 
the jets in microquasars are as relativistic as AGN 
jets provided the jets in microquasars are not confined 
and the derived opening angles are solely due to the 



relativistic effect llMiller- Jones. Fender, fc Naker 1 120061: 
iGopal-Krishna. Dhurde. fc Wiitall2004D . It is natiiral to 
associate (though disputed!) the radio loud and radio 
quiet dichotomy observed in AGN with jet-producing 
low hard stat e and non-jet-produ cing high soft state in 
microquasars l)Maccarone. Gallo. fc Fend er 2003). The 
PDS of luminous narrow line Seyfert 1 AGNs appears 
to be similar to the one observed during the high soft 
state in microquasars while the PDS of the low lu- 
minosity AGNs is similar to the one obs erved dur- 
ing the low hard state in microq uasars t|M'^Hardv et all 
l20f)l iMarkowitz fc Uttlev ll20(il . Shocks are related to 
the most energetic processes in the universe (Gamma 
ray bursts, large superluminal jets in microquasars and 
AGNs) and microquasars provide a great opportunity 
to study them under the best possible conditions within 
reasonable time. 

4. CONCLUSIONS 

We present for the first time clear evidences of direct 
connection between the accretion disk during plateaux 
and the following superluminal radio flares. We find that 
the rhaccr, as inferred from the X-ray fiux, is very high 
during the radio plateaux and hboi approaches I-iEdd- We 
suggest that such hot accretion disk during the radio 
plateaux always accompany with radiation-driven wind. 
The internal shock forms in the previously generated 
slowly moving wind from the accretion disk with /3 < 
0.01 as the fast moving discrete jet with P ^ 1 catches 
up and interacts with it. Both the components; slow 
moving wind and fast moving jet are related to the accre- 
tion disk during plateau state and the strength fc speed 
of these two should determine the power of the inter- 
nal shock. In GRS 1915-1-105 , it is a efficient internal 
shock scenario where not just velocity but total jet power 
has significantly increased. The superluminal flares have 
steep radio spectrum and are observed at large distances 
from the compact object (> 240 AU). The profile of the 
superluminal radio flares can be described as fast rise 
and exponential decay (FRED). 

The other radio flares observed in GRS 1915-1-105 in- 
clude the persistent radio emission during the plateaux 
(20-160 mJy), the radio oscillations fc discrete radio 
flares (5-150 mJy) and the preplateau radio flares 
(50-360 mJy). All of these are observed close to the 
compact object (< 200 AU) and have flat or inverted 
radio spectra. During the radio plateaux, both the 
radio flare and the accretion disk are in steady state 
(outflow and inflow in the accretion disk are in equi- 
librium) and a tight correlation is observed between 
the radio luminosity and the X-ray luminosity as 
discussed earlier. The IR & radio oscillations (5-150) 
are periodic ejections of adiabat ically expandi n g self 
absorbing synchrotron clouds IjMirabel et alJ Il998t 
llshwara-Chandra. Yadav. fc Pramesh Bad l2002tl . The 
IR flux suggests that the ejected matter decouples from 
the accretion disk at the time of the spike in the X-rays 
l)Eikenberrv et al.l Il998() . The preplateau radio flares 
(50-360 mJy) are discrete ejections which are closely 
spaced in time and hence produce overlapped radio flares. 
The preplateau flares have been modeled as adiabatically 
expanding self absorbing clouds ejected from the accre- 
tion disk which explain reasonably well all the available 
observed data of time delay for radio emission at lower 
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frequ ency Ijlshwara-Chandra. Yadav. fc Pramesh Raa 
I2002J) . It is widely believed that it is the coro- 
nal material and not the disk material which is 

ejected prior t o rad i o flares llRa u & Grciner 2003; 
^dawalc ct al. 2003: Fender. BcUoni. & Gallo 2004; 
Rothstein. Eikcnbcrrv. & Matthews 2005.). During 



oscillations/discrete flares and the prcplateau flares, the 
coronal mass ejections affect the hard X-ray component 
(Compton component) either partly or completely (state 
change) (outflow and inflow are not in equilibrium in 
this case unlike during the plateaux). 
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TABLE 1 

All selected superluminal radio flares and their properties (see text for details). X-ray Properties from RXTE 

PCA/HEXTE data during preceding radio plateau. 



Supcrluminal Rac 


io Flare Properties 




Associated Preceding Plateau Properties | 


M,JD/ 


Peak 


Rise 


Decay Time 


Radio 


GBI 


ASM 


Date of 


QPO 


Total 


Nh 


Date 


Flux 


Time 


constant 


Flux" 


Flux 


Flux 


RXTE 


Freq. 


X-ray 


(10^^ cm-2) 




(mJy) 


(day) 


(days) 


(Jy.day) 


(mjy) 


(c/s) 


Obser. 


(Hz) 


Flux'' 




50750/1997 Oct 30 


550 


<0.6 


3 16+"-"' 

■^■-'-"-0.09 


2.20* 


47.3 


35.9 


1997 Oct. 25 


1.88 


2.04 


-'-■^■^'-0.91 


50916/1998 Apr 13 


920 


<0.8 


, „2 + 0.36 


3.04 


91.0 


48.2 


1998 Apr. 11 


1.60 


2.42 


TA nr; + 0.S2 
14.95_o 55 


50933/1998 Apr 30 


580 


<0.9 


o no + 0-13 
d.9»_o 13 


3.02 


91.0 


44.9 


1998 Apr. 28 


1.41 


2.34 


,o r,-+0.78 
10.00_Q 3j^ 


50967/1998 Jun 03 


710 


0.3 


2 82+°-"'** 


2.15 


56.2 


36.5 


1998 May 31 


1.76 


2.14 


I2.99tlf, 


51204/1999 Jan 30 


340 


0.25 


, ,0 + 0.03 
-^•-^^-0.05 


0.36 


27.8 


33.8 


1999 Jan. 24 


2.55 


1.84 


"■i9;S:f5 


51337/1999 Jun 08 


490 


<0.7 


2 fi7+0-22 

^■°'-0.16 


1.58 


45.5 


35.2 


1999 Jun. 03 


1.77 


1.98 


12.43;?:- 


51535/1999 Dec 23 


510 


<0.8 


2 67+°-^^ 

^■°'-0.14 


1.75* 


51.3 


38.9 


1999 Dec. 21 


2.12 


2.23 


i3.6i;iii 


52105/2001 Jul 16"= 


210 


<0.7 


, 7- + 0-20 
'-■' '-0.20 


0.62* 


20.0 


30.0 


2001 Jul. 11 


2.40 


1.75 


io.6oiS:Io 



a: Integrated flux is obtained by fitting an exponential 
of Sxdccay time constant. Typical errors are ±0.05 
data or due to contribution by additional radio flare 

b: Integrated 3—150 keV X-ray flux in 10~ ergs cm 

c: VLBA radio data (see text for details) 



function to the flare proflle and integrating the function over duration 
(larger error (0.15—3.0) when marked "^" due to large gap in GBI radio 

) 



TABLE 2 
Best fit spectral parameters for. the X-ray observations prior to the radio flares shown in Table 1. 



Day of 




3 - 150 keV flux (in 10" 


-S -2 

ergs cm s 


"') kT.„ 


kTe 


r 


P 


RL 


x"(dof) 


RXTE obs. 


Ftotai 


Fdbfa 


Fatt 


Fpo. 


keV 


keV 






km 




1997 Oct. 


25 


2.04 


0.20 


0.75 


0.97 


0.50±:-'^^ 


5.50i:i^3 


7 7+O.S 

'■'-0.6 


9 47+0.30 
^■^' -0.03 


748.1 


1.40(74) 


1998 Apr. 


11 


2.42 


0.29 


0.93 


1.05 


0A6tlZ 


4.04«„;fo 


9 7+0.9 

T9-g 


2 07+0.09 
^■9'^-0.03 


1392.4 


0.83(94) 


1998 Apr. 


28 


2.34 


0.29 


1.04 


0.85 


o.5o;°:- 


4 i8+"-i^ 


a c+1-0 
^■8-o.4 


2 ^0+"-^^ 
^■OU-0.19 


831.4 


1.43(94) 


1998 May 


31 


2.14 


0.25 


0.85 


0.91 


o-5it.r„i 


4 96+°-^^ 


o r + l.O 
^■5-0.7 


2 (!n + 0.24 
^•69-0.06 


741.7 


0.93(94) 


1999 Jan. 


24 


1.84 


0.128 


0.60 


1.0 


0.57t°o:S 


7 00+°-^^ 

'■""-0.60 


4 2+1-9 


2 oq+0.14 
^••^^-0.09 


297.0 


1.15(93) 


1999 Jun. 


03 


1.98 


0.23 


0.73 


0.92 


tsi+o"i 
tSSA 


r no+0.80 


7 7+0.7 


2 ro + 0.23 


660.9 


1.16(88) 


1999 Dec. 


21 


2.23 


0.25 


0.71 


1.16 


0.48t°o:S? 


- „c:+0.76 
"■^°-0.39 


8.3+^-^ 


o OO + 0.03 


981.3 


0.81(88) 


2001 Jul. 


11 


1.75 


0.13 


0.50 


0.95 


0.54+°:"? 


6.99l°o:f9 


7.it';:=5 


2 on + 0.03 
^■•59-0.07 


354.3 


0.98(87) 



'': Rin derived assuming distance D — 12.5 kpc and inclination 
associated best fit values of N/f is shown in Table 1 



angle, 



: 70° 



